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Abstract 

We study the oblate-prolate shape mixing in the low-lying states of proton- 
rich Kr isotopes using the five- dimensional quadrupole collective Hamilto- 
nian. The collective Hamiltonian is derived microscopically by means of the 
CHFB (constrained Hartree-Fock-Bogoliubov) + Local QRPA (quasiparticle 
random phase approximation) method, which we have developed recently on 
the basis of the adiabatic self-consistent collective coordinate method. The 
results of the numerical calculation show the importance of large-amplitude 
collective vibrations in the triaxial shape degree of freedom and rotational 
effects on the oblate-prolate shape mixing dynamics in the low-lying states 
of these isotopes. 
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1. Introduction 

Atomic nuclei exhibit different shapes depending on the numbers of pro- 
ton and neutron, the excitation energies, or angular momentum. Shape co- 
existence phenomena, in which an excited band with a shape different from 
the shape in the ground band exists close to the ground band, are widely 
observed all over the nuclear chart. From the mean-field viewpoint, shape 
coexistence indicates that there are two equilibrium points in the mean field, 
and appearance of a low-lying excited 0^ state is one of typical indications 
of the shape coexistence. 

In proton-rich Kr isotopes, it has been conjectured that the oblate and 
prolate shapes coexist in their low-lying states, since the ground bands quite 
different from regular rotational spectra [1, 2] and the low-lying excited O"*" 
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states [3, 4] were measured. Since then, much experimental data have been 
accumulated for proton-rich Kr isotopes [5, 6, 7, 8] and they support the 
interpretation as the oblate-prolate shape coexistence. The spectroscopic 
quadrupole moments [6] suggest the prolate ground state in ^^Kr and ^^Kr, 
while the ground state of ^^Kr is assumed to be oblate from the properties of 
the E2 transition probabilities [7] . The systematics of the electric monopole 
transition strengths [3, 4, 9] and the excitation energies of the first excited 
0"*" states are consistent with the interpretation that a shape transition from 
the oblate ground state in ''^Kr to the prolate ground states in ''^Kr and ^^Kr 
takes place. 

The description of these isotopes can be a touchstone for nuclear structure 
models, and many attempts were done from several theoretical approaches 
[10, 11, 12, 13, 14, 15, 16]. Recently Girod et al. [17] have reproduced the 
shape transition from the oblate ground state in ''^Kr to the prolate in ^^Kr 
with the HFB (Hartree-Fock-Bogoliubov)-based GCM (generator coordinate 
method) + GOA (Gaussian overlap approximation) calculation. Making a 
comparison with the result of the axial GCM calculation done by Bender et 
al. [15], they have pointed out that the triaxial shape plays an essential role 
in the shape coexistence and the shape transition in the light Kr isotopes. 
The importance of large-amplitude collective vibrations in the triaxial shape 
degree of freedom has also been shown in the previous works for ^^Kr [18] 
and neighboring ^^"''^Se [19] on the basis of the (1+3) dimensional calculation 
using the adiabatic self-consistent collective coordinate (ASCC) method [20, 
21]. 

In this paper, we study the low-lying states of proton-rich Kr isotopes us- 
ing a new method we have developed recently [22] . In this method, we deter- 
mine the five-dimensional (5D) quadrupole collective Hamiltonian [23, 24, 25, 
26] microscopically on the basis of the ASCC method. Local normal modes 
on top of each constrained HFB (CHFB) state at each point on the (/3,7) 
plane is the main concept. After solving the CHFB equations imposing the 
constraints on the deformation and the particle numbers, we solve the local 
QRPA (LQRPA) equations, which is an extension of the usual QRPA (quasi- 
particle random phase approximation) to non-HFB-equilibrium points, on 
top of the CHFB states. Therefore, we call this method the CHFB+LQRPA 
method. Using this method, we derive the seven quantities in the collective 
Hamiltonian: the collective potential, three vibrational inertial masses, and 
three rotational moments of inertia. 

One of the advantages of the CHFB+LQRPA method is that the vibra- 
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tional and rotational masses (inertial functions) determined with this method 
contain the contributions from the time-odd components of the mean-field 
unlike the widely-used Inglis-Belyaev (IB) cranking masses [27, 28]. It is 
well known that the ignorance of the time-odd components in the IB masses 
breaks the self-consistency of the theory [29, 30]. Nevertheless, even in recent 
microscopic studies by means of the 5D quadrupole collective Hamiltonian 
[17, 31, 32, 33, 34], the IB cranking masses are still used at least for vibra- 
tional inertial functions. In this study, we use the pairing-plus-quadrupole 
(P-l-Q) force model [35, 36] including the quadrupole-pairing force and take 
into account the contributions from the time-odd components of the mean 
field to the vibrational and rotational masses. Inclusion of the quadrupole 
pairing force is essential because it is the only term which gives the time-odd 
components of the mean field [37]. 

This paper is organized as follows. In Section 2, we summarize the proce- 
dure of deriving the 5D quadrupole collective Hamiltonian by means of the 
CHFB+LQRPA method. In Section 3, we calculate the vibrational and ro- 
tational masses. By solving the collective Schrodinger equation, we calculate 
excitation spectra, B{E2), spectroscopic quadrupole moments and monopole 
transition matrix elements for the low-lying states in '^2,74,76j<^j, ^^^^ discuss 
the properties of the oblate-prolate shape coexistence/mixing in these nuclei. 
Conclusion is given in Section 4. 

2. Microscopic derivation of the 5D quadrupole collective Hamil- 
tonian 

In this section, we summarize the procedure of microscopically deriv- 
ing the 5D quadrupole collective Hamiltonian and the collective Schrodinger 
equation by means of the CHFB+LQRPA method. 

2.1. 5D quadrupole Hamiltonian and collective Schrodinger equation 

The 5D quadrupole collective Hamiltonian is written in terms of the mag- 
nitude f3 and the degree of triaxiality 7 of quadrupole deformation, and their 
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time derivatives, (3 and 7, as 

"HcoU 



Tvib+Trot + V^(/3,7), 

Tvib = ^/^/3/3(/3,7)/3' + I)/37(/3,7)/37 



^77(/3,7)7', 
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T™t = ^$^Jfc(/3,7) 



2 



(1) 

(2) 
(3) 



k=i 



where Tyib, T^ot and represent the vibrational, rotational and collective 
potential energies, respectively. The 5D quadrupole collective Hamiltonian 
has seven quantities to be determined: three vibrational inertial masses, three 
rotational moments of inertia, and the collective potential. The moments of 
inertia can be parametrized as 



7) = A(3^Dk{(3, 7) sin2 7, {k = 1, 2, 3) 



(4) 



with 7fc = 7 — (27rA;/3). The vibrational and rotational masses are, in general, 
functions of the deformation parameters. How to determine these inertial 
masses and the collective potential is discussed in the following subsection. 

As the collective Hamiltonian (1) is classical, we quantize it according to 
the Pauli prescription to obtain the collective Schrodinger equation 



where 



{TVib + Trot + ^}^a7M(/3, 7, ^) = ^a/^a/Af(/3, 7, ^) 
-1 f 
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— sin d,'-fDp^dp I + 
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- rot 



2Jfc' 



Here, R{(3,'j) and W{f3,'y) are defined as 

i?(/3,7) =Z}i(/3,7)D2(/3,7)^3(/3,7), 
W{(3,j)={D0p{(3,-f)D,,i(3,^) - [D;3,(/3,7)]'}/3" 



(5) 

(6) 
(7) 

(8) 



(9) 
(10) 
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The collective wave function \I/q,/m(/5, 7, ^) is specified by the total angu- 
lar momentum /, its projection onto the 2;-axis of the laboratory frame M, 
and a distinguishing the states with the same I and M. As the collective 
potential is independent of the Euler angles Q, the collective wave function 
is written in the form: 



^o.IM{/3,l,n)= ^aIK{P.l){WMK), (11) 



_fsr=even 



where 



{n\iMK) = ^ —^±l^^[vi,j,{n) + {-YviMm- (12) 

Here, V^^j^ is Wi gner's rotation matrix and K is the projection of the angular 
momentum onto the 2;-axis in the body-fixed frame. The summation over K 
is taken from to / for even I and from 2 to J — 1 for odd /. 

The vibrational wave functions in the body-fixed frame, $a/x(/3,7), are 
normalized as 



^ ci/3ci7l$a/(/3,7)r|G'(/3,7)|^ = l, (13) 
where 

|$„,(/9,7)r= E \^o.iK{f3,7)\', (14) 

K=even 

and the volume element G(/3,7) is given by 

G{f3, 7) = 4{3^W{f3, ^)R{f3, 7) sin^ 87. (15) 

The symmetries and boundary conditions of the collective Hamiltonian and 
wave function are discussed in Ref. [25]. 

2.2. Constrained HFB + Local QRPA method 

In this subsection, we summarize the method of determining the inertial 
masses and collective potential. This method can be regarded as an approx- 
imation of the two-dimensional (2D) version of the ASCC method. In this 
method, we solve the local QRPA equations at each point on the (/3, 7) plane 
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on top the CHFB state, and therefore we call it the CHFB+LQRPA method. 
(See Ref. [22] for details.) 

First, we solve the CHFB equation given by 



6 (0(/3, 7) I ^chfb(/3, 7) 7)) = 0, (16) 



with four constraints 

(0(/3,7)|iVW |0(/3,7)) = (r = n,p) (18) 

(0(/3,7)|^Sl0(/3,7)) = ^$l\ (m = 0,2) (19) 

where 1)2^ denotes Hermitian quadrupole operators, D20 and (Z)22 + -D2_2)/2 
(for m = and 2, respectively). We define the quadrupole deformation vari- 
ables (/3, 7) in terms of the expectation values of the quadrupole operators: 

/3cos7 = r/D^+\ (20) 
-^/3sin7 = r/D£\ (21) 

where 77 is a scaling factor (to be discussed in Section 2-4). 

Second, we solve the LQRPA equations on top of the CHFB states ob- 
tained above. 



5(0(/3,7)|[^chfb(/3,7),QX/3,7)] 
-iA(/3,7)l0(/3,7))=O, (22) 

<5(0(/3,7)|[^chfb(/3,7)4A(/3,7)] 
- a(/3, l)Q\f3, 7) m, 7)) = 0, (^ = 1, 2). (23) 

Here the infinitesimal generators, Q*(/3,7) and Pj(/3,7), are local operators 
defined, for a given set of quadrupole deformation variables (/3,7), with re- 
spect to the CHFB state |0(/3,7)). The quantity Ci(/3,7) is related to the 
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eigenfrequency 7) of the local normal mode through uf{f3, 7) = Ci{f3, 7). 
Note that these equations are valid also for regions with negative curvature 
(Ci(/3,7) < 0) where Ui{P,'-f) takes an imaginary value. For selecting two 
collective modes from among many LQRPA modes, we use the criterion for- 
mulated in Ref. [22]. 

Using the time derivatives of D^, the collective vibrational energy (2) 
is written in a form of 

Tvib = ^MooiDiPf + Mo2DiVD^2^^ + ^Masp^? )^ (24) 

In the 2D ASCC method, the vibrational part of the collective Hamiltonian 
is given by 

i 

under the condition that there exists a collective coordinate system {qi,q2) 
where the vibrational masses and stiffness tensors can be diagonalized glob- 
ally. Assuming one-to-one correspondence between (gi,g2) and (-D2o'', -D22 
the vibrational masses in Eq. (24) are written as 

i=l,2 '^-'^2m '~^^2m' 

In this approximate version of the 2D ASCC method, we equate the momen- 
tum operator in the ASCC method with the one obtained by the LQRPA 
equations, and then the partial derivatives are easily evaluated as 

= [B£',iA(A7)l m,!)), (27) 

^ 

^{4>{i3,7)\ ib'iK-m,',)] m,-/)), (28) 

without need of numerical derivatives. 

The rotational moments of inertia are calculated by solving the LQRPA 
equation [22] for rotation on top of each CHFB state: 
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<5(0(/3,7)| [i^cHFB,^fc(/3,7)] - -(Jfe)-^4|0(/3,7)) = 0, (29) 
(0(/3,7)| [^fc(/3,7),4,] |0(/3,7)) = (30) 

where \l/fc(/3,7) and represent the rotational angle and the angular mo- 
mentum operators with respect to the principal axes associated with the 
CHFB state |0(/3,7)). This is an extension of the Thouless-Valatin equa- 
tion [38] for the HFB equilibrium state to non-equilibrium CHFB states. We 
call i7fc(/3,7) and Dfc(/3,7) determined by the above equations and Eq. (4) 
'LQRPA moments of inertia' and 'LQRPA rotational masses', respectively. 

Last, we solve the collective Schrodinger equation with the 5D quadrupole 
collective Hamiltonian whose collective potential and inertial masses are cal- 
culated as above. 

2.3. Calculation of electric quadrupole (E2) transitions and moments 
The E2 operator in the intrinsic frame is defined as 

= E ^iJ^S, (31) 

r=n,p 

where D^2, quadrupole operator for protons or neutrons, and thus 

X]r-^2m — D2m- The E2 Operator in the laboratory frame is given by 

^r^ = E^w(^)^S'' (32) 

m' 

with Wigner's V functions. The reduced E2 transition probability B[E2) 
and the spectroscopic quadrupole moment Q are given by 

B{E2-aI ^ a'l') = {21 + ly^ \{aI\\D'^''^\\a' (33) 
Q{al) = (a, /, M = /| dJ^'^ |a, /, M = /) (34) 

The reduced matrix element in Eq. (33) is easily evaluated through the 
relation, 
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Substituting Eq. (11) into M), we obtain 



V(2/+l)(2/' + l)(-/5^|(/ ^ ^' ^''^ --'^^""'^ 



+ 



/ 2 1' 

K-2 2 A' 

12 r 

K 2 -K-2 



^a,I,K+2 \ -^2+^^ \^a',I',K) 
^2+ \^a',I',K+2) 



(36) 



with Dif^ = (D^^ + dSV2. 

The intrinsic matrix elements appearing in the above formula are evalu- 
ated as 



{^aj,K\Di^+\^a',i',K') = I c?/3d7|G|^$U(/3,7)^mT(/3,7)$a'/'K'(/3,7) 



(E2), 



where 



(S2) 



(37) 
(38) 



We also calculate the electric monopole transition matrix elements. It is 
evaluated, to the lowest order in /3, by [39]. 



p{E0;z^f) = l-{^(3'\f)6j^,^ 

4 71 

(39) 

2.4- Details of numerical calculation 

In this work, we adopt a version of the P+Q interaction which includes 
the quadrupole pairing interaction as well as the monopole pairing interac- 
tion. We take two major shells as the active model space for protons and 
neutrons. The interaction parameters are determined as follows: for ''^Kr, we 
use the same values of the monopole pairing interaction strength G^^^ and the 
quadrupole particle-hole interaction strength x as used in the previous ASCC 
work [18], which were determined such that the magnitude of the quadrupole 
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deformation /3 and monopole pairing gaps at the oblate and prolate HFB 
minima obtained in the Skyrme-HFB calculation [14] were approximately re- 
produced. These values are scaled for ^^Kr and ^^Kr assuming a simple mass 
number dependence [36]: G^^^ ~ A~'^ and x' = X^^ ~ A~^. We determine 
the quadrupole-pairing strength following the Sakamoto-Kishimoto prescrip- 
tion [40]. As in the conventional treatment of the P+Q model, we omit the 
Fock term. Therefore, we use the abbreviation HB instead of HFB in the 
following. The CHB + LQRPA equations are solved at each mesh point on 
the (/3,7) plane. This two-dimensional mesh consists of the 60 x 60 points 
specified by 



with which we can perform a parallel computation easily for each mesh point. 
The effective charges are given by e^^^ = 6epo\ for neutrons and = e + 
Scpoi for protons. In this paper, we use the polarization charge 6epo\ = 
0.834e, which is determined such that the calculated value oi B(E2; 2f or) 
reproduces the experimental data for ''^Kr [6]. 

3. Results of calculation and discussion 

In this section, we present the result of the numerical calculation for ''^Kr, 
^■^Kr, and ^^Kr and discuss properties of low-lying collective modes of exci- 
tation in these nuclei. Figure 1 shows the collective potential calculated for 
72-76j<^j^ All of the potential have two local minima: one is oblate and the 
other is prolate. The spherical shape is a local maximum in any case. In 
^^Kr, the absolute minimum is oblate and the prolate minimum is approxi- 
mately 600 keV higher. In contrast, the prolate minimum is lower than the 
oblate one in ^^Kr, which suggest there may occur a shape transition from 
oblate to prolate between these isotopes. The experimental data [6] of the 
spectroscopic quadrupole moments for ^^Kr and ^^Kr indicate that they are 
prolate in the ground band. One may expect that the absolute minimum is 
oblate in ^^Kr and it becomes prolate for A > 74. However, it is not as simple 
as one might expect: the oblate minimum becomes the absolute minimum 
again in ^^Kr. This does not necessarily mean that our calculation fails in the 
reproduction of the ground band shape in ''^Kr. As we shall discuss below, 
we have to take into account dynamical effects beyond the mean field. 



= (i- 0.5) X 0.01, (i = l,---60), 
7, = (j - 0.5) X 1°, (j = l,...60). 



(40) 
(41) 
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We show in Fig. 2 the /3 — 7 dependence of the vibrational masses for Kr. 
One can see their deviation from a constant. Although they have a compli- 
cated structure at large deformation compared to that at small deformation, 
this region hardly contributes to the calculated results as the potential is 
high. We also depict the ratios of the LQRPA vibrational masses to the IB 
vibrational masses in Fig. 3. 

We plot the LQRPA moments of inertia for ^^Kr in Fig. 4. We can see 
deviation from the irrotational moments of inertia. The ratios of the LQRPA 
moments of inertia to the Inglis-Belyaev ones are plotted in Fig. 5. All the 
ratios shown in Figs. 3 and 5 are larger than 1 on the entire {(3, 7) plane and 
depends on /3 and 7. This indicates not only that the time-odd components 
of the mean field (which enhance the inertial masses) should be taken into 
account but also that a simple remedy used in Refs. [31, 32, 33], i. e., a 
simple multiplication of a constant factor to the IB cranking moments of 
inertia, is insufficient. 

3.1. ^^Kr 

We show the excitation spectra and B{E2) values calculated for '^^Kr in 
Fig. 6 together with the experimental data. The result of calculation for spec- 
troscopic quadrupole moments are shown in Fig. 13 together with those for 
74,76j^j,^ In Fig. 7 the collective wave function squared, /3^J2k \^aiK{/3,l)\'^, 
are plotted with a factor multiplied. The (3'^ factor carries the dominant 
(3 dependence of the volume element G^/^(/3,7). Note that, if the whole vol- 
ume element was multiplied, all the state would look triaxial because of the 
factor sin 37, which vanishes on the oblate and prolate lines. In Fig. 6, the 
result obtained using the IB cranking masses instead of the LQRPA collective 
masses are also shown for comparison. 

As expected from Figs. 3 and 5, the excitation energies obtained with 
the IB cranking masses are higher than those obtained with the LQRPA col- 
lective masses. Moreover, the excitation energies obtained with the LQRPA 
collective masses are in better agreement with the experimental data except 
for the 2f state. The observed excitation energy of the O2 state locates close 
to the 2f state. It is seen that the 0^ energy obtained with the LQRPA 
masses is in much better agreement with the experimental data than that 
with the IB cranking masses. 

The collective wave function of the ground state has a peak around the 
oblate potential minimum. It has a tail to the prolate region, however. As 
angular momentum increases, the localization of the collective wave function 
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develops. The collective wave function of the 0^ state consists of two com- 
ponent: one is a sharp peak on the oblate side and the other is a component 
spreading around the prolate region somewhat broadly. It has a node in the 
P direction on the oblate side but the 7-vibrational component is strongly 
mixed. In contrast to the 0^ state, the effects of the f3 vibrational excitation 
are weak for the other yrare states. The vibrational wave function has two 
peaks at / = 2. The oblate peak shrinks rapidly but the prolate one grows 
with angular momentum due to the orthogonality to the yrast states. These 
behaviors of the vibrational wave function are consistent with the results 
of the (1+3) dimensional ASCC calculation [18, 19]. There, it is suggested 
that the rotational effect may assist the localization of the collective wave 
function. 

The E2 transition strengths clearly indicate that the shape-coexistence- 
like character becomes stronger with increasing angular momentum: the in- 
terband transitions between the initial and final states having equal angular 
momentum become weaker and weaker. 

The signs of the spectroscopic quadrupole moments shown in Fig. 13(a) 
are positive for the yrast states indicating their oblate-like character, while 
they are negative for the yrare states indicating their prolate-like character. 
Their magnitudes increase as angular momentum increases, which reflects 
the growth of localization of the vibrational wave function. 

3.2. '^^Kr 

The excitation spectrum and B{E2) values calculated for ''^Kr are pre- 
sented in Fig. 8 together with the experimental data, while the collective 
wave function squared are displayed in Fig. 9. It is seen that the collective 
wave function of the ground state spreads over the entire 7 region along 
the potential valley. However, it localizes more and more on the prolate 
side with increasing angular momentum. This behavior results from the /3- 
7 dependence of the rotational moments of inertia: one can clearly see the 
oblate-prolate asymmetry for /3 > 0.4 in the moment of inertia J7i displayed 
in Fig. 10. Although the 0^ state has somewhat a /3- vibrational component, 
the yrare collective wave function for / > 2 has almost one-dimensional struc- 
ture in the 7 direction with a constant (3. While the collective wave function 
has a two-peak structure, the prolate peak shrinks with increasing angular 
momentum due to the orthogonality to the yrast state. 

The third and fourth states for each angular momentum can be regarded 
as admixtures of /3-vibrational and 7-vibrational components. On one hand. 
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the 0^ and 2j states have a node in the (3 direction on the oblate side. On 
the other hand, O4 , 23 , 4jj" and have a node on the prolate side. 

It is clearly seen in Fig. 8 that the result of calculation reproduces well 
the significant enhancement of i?(£'2;0j' — 2f), i?(_E2; 2^ — )■ 0^) and 
— )■ 2^), found in experiments, as well as the large B(^E2) values 
within the ground band. 

The spectroscopic quadrupole moments for ^^Kr, shown in Fig. 13(b), are 
also in excellent agreement with the experimental data, aside from a minor 
disagreement for the 2^ state. In particular, the characteristic features, such 
as their signs and the increasing tendency of their magnitudes with angular 
momentum in the ground band, are well reproduced. 

3.3. '^^Kr 

The excitation spectrum and B{E2) values calculated for ^^Kr are pre- 
sented in Fig. 11 together with the experimental data, while the collective 
wave function squared are displayed in Fig. 12. The behavior of the collec- 
tive wave function for ^^Kr is qualitatively similar to that for ''^Kr, except 
that the development of the localization of the yrare wave function is weaker 
and it has a more remarkable two-peak structure. The mechanism of ap- 
pearance of the two-peak structure is discussed from a general viewpoint of 
oblate-prolate symmetry breaking in Ref. [41]. 

The difference to be most noted between ^^Kr and ^^Kr is the location 
of the minimum of the collective potential. One can see the importance of 
the rotational effect from Fig. 12. Although the potential has the oblate 
absolute minimum in ^^Kr, the vibrational wave functions of the 2f, 4f 
and 6f states still localize on the prolate side. It implies that the potential 
and rotational effects compete and the rotational effect dominates over the 
potential one with increasing angular momentum. This calculation serves as 
a good example suggesting that the nuclear shape cannot be determined only 
by the location of the potential minimum. 

The spectroscopic quadrupole moments, shown in Fig. 13(c), qualitatively 
agree with the experimental data. The sign and increase with angular mo- 
mentum are well reproduced for the ground band. Their absolute magnitudes 
are, however, too small compared with the experimental data. This seems 
to indicate that the calculated magnitudes of the quadrupole deformation 
and/or the localization of the collective wave functions in the (/3,7) space 
are insufficient. It is interesting to note that the calculated quadrupole mo- 
ment for the 2^ state better agrees with the experimental data for the 2^ 
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state. In this connection, we notice that the E2 transition properties of the 
observed 2^ state in ™Kr are different from those in '''^Kr: while it decays to 
the 2^ and 0^ states with almost equal E2 strengths in ''^Kr, it decays al- 
most totally to the 0^ state in ''^Kr. This point was emphasized by Clement 
et al. [6]. This difference is qualitatively reproduced in our calculation if we 
assume that, in ^^Kr, the band consisting of the experimental 0^, 2^, 4^ and 
6^ states corresponds to the band comprised of the calculated 0^, 23 , 43 and 
63 states. 

3.4- Discussion 

In addition to the E2 transition properties and the quadrupole mo- 
ments of the low-lying states, we have also calculated EO transition strengths 
p^(-EO). The result of calculation is compared with the experimental data 
in Fig. 14. We see that the calculated result well reproduces the experimen- 
tal data both qualitatively and quantitatively. The p'^{EO) takes a maximal 
value a.t A = 74, which reflects the shape transition. These EO transition 
strengths are evaluated using Eq. (39) where the phenomenological Bohr rela- 
tion [23, 39] between the £"0 matrix elements and the quadrupole deformation 
is assumed. In this sense, this way of calculating p{EO) is semi-microscopic, in 
contrast to the calculation of B{E2) and spectroscopic quadrupole moments 
which are carried out in a fully microscopic way using the expectation values 
of the quadrupole operators at each point on the 7) deformation space, as 
explained in Section 2.3. For a fully microscopic calculation of p{EO) without 
assuming Bohr's relation, we need to use a more realistic effective interaction 
than Baranger-Kumar's version of the P+Q force model [36]. This task is left 
for future works. Note that p{EO) vanishes in Baranger-Kumar's P+Q force 
model, because this force is designed such that the radial matrix elements of 
the monopole operators take the same value for all single-particle states. 

Let us add a few remarks concerning the results of calculation for E2 
transition properties shown in Figs. 8 and 11. Although we have succeeded 
in reproducing a number of features observed in experiments, some aspects 
in experimental data remains unexplained. In particular, the large value of 
B{E2] 0^ — )■ 2^) in ^^Kr is not reproduced in the calculation. In this con- 
nection, it should be emphasized that the properties of the excited 0"'' states 
are quite sensitive to the mixing of the (3 vibrations and the large-amplitude 
7 vibrations (large-amplitude fluctuations in the triaxial shape degree of 
freedom). In this study, we have not adjusted the interaction parameters by 
fitting the theoretical result to the experiment, and the interaction we have 
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employed itself is rather simple. Therefore, it remains for future to examine 
whether or not these disagreements with experiments are improved when a 
better interaction is used. 

As we mentioned in Introduction, the low-lying states of light Kr isotopes 
have recently been studied also by Girod et al. [17] by means of the CHFB- 
based GCM (GOA) method. Consistent with our results, their calculation 
also demonstrates that the large-amplitude fluctuation in the triaxial shape 
degree of freedom plays an important role in the description of the shape co- 
existence and the shape transition in the light Kr isotopes. They have used 
the Gogny DIS force, which is more realistic than the P-l-Q force. However, 
their inertial masses are derived with the cranking formula, which ignores 
the contribution from the time-odd mean field. We would like to empha- 
size that there is no problem to use a more realistic interaction within the 
framework of the CHFB+LQRPA method, although numerical calculation 
becomes much more demanding. The reason why we have used the schematic 
P+Q interaction in this paper is simply because this is the first application 
of the CHFB-I-LQRPA method to real nuclear collective phenomena. We 
shall employ a more realistic interaction like Skyrme interactions in future 
works. It will also become possible to use the CHFB+LQRPA method in 
conjunction with modern nuclear energy density functionals. 

4. Conclusion 

We have studied the oblate-prolate shape coexistence/mixing in the low- 
lying states of '^2,74,76j^j, using the 5D quadrupole Hamiltonian determined 
microscopically by means of the CHFB+LQRPA method, which is based on 
the 2D ASCC method. Our results indicates a shape transition from oblate 
in the ground state of ''^Kr to prolate in ''^'''^Kr, which is consistent with the 
experimental data. We have shown that the basic features of the low-lying 
spectra in these nuclei are determined by the interplay of the large-amplitude 
vibrations in the triaxial shape degree of freedom, the /3-vibrational excita- 
tions and the rotational motions. We have furthermore shown that the ro- 
tational motion plays an important role for the growth of localization of the 
vibrational wave functions in the 7) deformation space. 
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discussions. K. S. and N. H. are supported by the Junior Research Associate 
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(a) 72Kr (b) ^^Kr (c) ^^Kr 



Figure 1: Collective potential energy surfaces V{l3,j) for "^2,74,76]^^.^ r^j^^ regions higher 
than 5 MeV (measured from the HFB minima) are colored rosy-brown. 
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Figure 2: Vibrational inertial masses, _D/3^(/3,7), Z?^^(/3,7)//3 and D^j{f3,j)/f3'^, in unit 
of MeV"^ calculated for '^^Kr. 
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(a) Ji(/3,7) (b) J2(/3,7) (c) MP.l) 

Figure 4: LQRPA rotational moments of inertia in unit of MeV~^, calculated for ^^Kr. 




(a) Ji/J['''^ (b) J./J^'''^ (c) Ja/J-i^^) 



Figure 5: Ratios of the LQRPA rotational moments of inertia to the IB cranking moments 
of inertia, calculated for ^^Kr. 
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Figure 6: Excitation spectra and B{E2) values calculated for ^^Kr by means of the 
CHB+LQRPA method (denoted LQRPA) and experimental data [4, 7, 8]. For com- 
parison, results calculated using the IB cranking masses (denoted Inglis-Belyaev) are also 
shown. Only i?(i?2)'s larger than 1 Wcisskopf unit are shown in units of e^fm*. 





Figure 7: Vibrational wave functions squared, /3^|$q,7(/3, 7)^, for "^^Kr. 
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Figure 9: Vibrational wave functions squared, /3'*|$q,7(/3, 7)^, for ^''Kr. 
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Figure 10: The LQRPA moment of inertia J\ for rotation about the x-axis, calculated for 
74Kr. 
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Figure 12: Vibrational wave functions squared, /3 |$q/(/3,7)| , for Kr. 
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Figure 13: Spectroscopic quadrupole moments in unit of e fm^ of the first (square), second 
(circle) and third (triangle) states for each angular momentum in 'i'2.74,76j^j._ Calculated 
values are shown by open symbols, while experimental data [6] are indicated by filled 
symbols. 
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Figure 14: Electric monopolc transition matrix element p^{EQ) calculated using Eq. (39) 
in comparison with the experimental data [3, 4] and the result of the HFB-based GCM 
(GOA) calculation [17]. 
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